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Abstract. The Travelling Salesman Problem (TSP) is a well known 
and challenging combinatorial optimisation problem. Its computational 

intractability has attracted a number of heuristic approaches to generate 
satisfactory, if not optimal, candidate solutions. Some methods take their 
inspiration from natural systems, extracting the salient features of such 
systems for use in classical computer algorithms. In this paper we demon- 
strate a simple unconventional computation method to approximate the 
Euclidean TSP using a virtual material approach. The morphological 
adaptation behaviour of the material emerges from the low-level inter- 
actions of a population of particles moving within a diffusive lattice. A 
'blob' of this material is placed over a set of data points projected into 
the lattice, representing TSP city locations, and the blob is reduced in 
size over time. As the blob shrinks it morphologically adapts to the con- 
figuration of the cities. The shrinkage process automatically stops when 
the blob no longer completely covers all cities. By manually tracing the 
perimeter of the blob a path between cities is elicited corresponding to 
a TSP tour. Over 6 runs on 20 randomly generated datasets of 20 cities 
this simple and unguided method found tours with a mean best tour 
length of 1.04, mean average tour length of 1.07 and mean worst tour 
length of 1.09 when expressed as a fraction of the minimal tour com- 
puted by an exact TSP solver. We examine the insertion mechanism by 
which the blob constructs a tour, note some properties and limitations 
of its performance, and discuss the relationship between the blob TSP 
and proximity graphs which group points on the plane. The method is 
notable for its simplicity and the spatially represented mechanical mode 
of its operation. We discuss similarities between this method and previ- 
ously suggested models of human performance on the TSP and suggest 
possibilities for further improvement. 
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1 Introduction 



The Travelling Salesman Problem (TSP) is a combinatorial optimisation prob- 
lem well studied in computer science, operations research and mathematics. In 
the most famous variant of the problem a hypothetical salesman has to visit a 
number of cities, visiting each city only once, before ending the journey at the 
original starting city. The shortest path, or tour, of cities, amongst all possible 
tours is the solution to the problem. The problem is of particular interest since 
the number of candidate solutions increases greatly as n, the number of cities, 
increases. The number of possible tours can be stated as (n — l)!/2 which, for 
large numbers of n, renders assessment of every possible candidate tour compu- 
tationally intractable. Besides being of theoretical interest, efficient solutions to 
the TSP have practical applications such as in vehicle routing, tool path length 
minimisation, and efScient warehouse storage and retrieval. 

The intractable nature of the TSP has led to the development of a number of 
heuristic approaches which can produce very short — but not guaranteed mini- 
mal — tours. A number of heuristic approaches are inspired by mechanisms seen 
in natural and biological systems. These methods attempt to efficiently traverse 
the candidate search space whilst avoiding only locally minimal solutions and 
include neural network approaches (most famously in [2S]), evolutionary algo- 
rithms [35], simulated annealing methods [23] i the elastic network approaches 
prompted in [16| . ant colony optimisation [13] . living [3] and virtual [29| slime 
mould based approaches, and bumblebee foraging [55] . 

Human performance on the TSP has also been studied in both naive and 
tutored subjects (see, for example, [35]). This is of particular interest because, 
unlike many nature inspired approaches, the human computation of TSP is by 
an individual and not based on population methods which evaluate a number 
of candidate solutions. Human performance on the TSP is also, for a limited 
number of cities at least, comparable in performance with heuristic approaches 
[21] , [M] . Although there are a number of competing theories as to how exactly 
humans approximate the TSP [37] , [H] , [33] , discovery of the methods employed 
may be useful as an insight into the mechanisms underlying complex percep- 
tual and cognitive processes and potentially as an aid for the development of 
computational algorithms. 

In this paper we adopt a material-based, minimum complexity approach. We 
show how a spatially represented non-classical, or unconventional, computational 
mechanism can be used to approximate the TSP. Taking inspiration from the 
non-neural, material-based computational behaviour of slime mould, we employ 
a sheet, or 'blob' of virtual material which is placed over a spatial map of cities. 
By shrinking this blob over time, it conforms and adapts to the arrangement of 
cities and a tour of the TSP is formed. We give an overview of the inspiration 
for the method in Section [2] The shrinking blob method is described in Section 
[3] Examples of the performance of the method compared to exact solutions 
generated by a TSP solver are given in Section [4] along with an analysis of the 
underlying mechanism and factors affecting the performance of the approach. 
We conclude in Section [5] by summarising the approach and its contribution in 



terms of simplicity. We examine similarities between the underlying mechanism 
of the shrinking blob method and proposed models of TSP tour perception and 
construction in studies of human performance on the TSP. We suggest further 
research aimed at improving the method. 

2 Slime Mould Inspired Computation of the TSP 

The giant single-celled amoeboid organism, true slime mould Physarum poly- 
cephalum, has recently been of interest as a candidate organism for the study of 
non-neural distributed computation. In the vegetative Plasmodium stage of its 
complex life cycle the organism forages towards, engulfs and consumes micro- 
organisms growing on vegetative matter. When presented with a spatial con- 
figuration of nutrient sources the Plasmodium forms a network of protoplasmic 
tubes connecting the nutrients. This is achieved without recourse to any spe- 
cialised neural tissue. The organism dynamically adapts its morphology to form 
efficient paths (in terms of a trade-off between overall distance and resilience to 
random damage) between the food sources |39) . [38| . |40| . 

Research into the computational abilities of Physarum was prompted by Nak- 
agaki, Yamada and Toth, who reported the ability of the Physarum Plasmodium 
to solve a simple maze problem |41) . It has since been demonstrated that the 
Plasmodium successfully approximates spatial representations of various graph 
problems |39I45I2I2"7| , combinatorial optimisation problems |7l8j , construction of 
logic gates and adding circuits |48I30I3) . and spatially represented logical ma- 
chines |l|5j . 

Although Physarum has been previously been used in the approximation of 
TSP [S] , this was achieved by an indirect encoding of the problem representation 
to enable it to be presented to a confined plasmodium in a controlled environ- 
ment. In the work by Aono et. al. it was shown that the morphology of the 
Plasmodium confined in a stellate chamber could be dynamically controlled by 
light irradiation of its boundary. When coupled to an elegant feedback mech- 
anism using an analysis method (to assess the presence of plasmodium at the 
extremities of the chamber), combined with Hopfield-Tank type neural network 
rules |25j . the plasmodium was used to generate candidate solutions to simple 
instances of the TSP |7|8j .In its natural propagative state, however, Physarum 
does not approximate area representations of a set of points, including the Con- 
vex Hull, Concave Hull [4 and the TSP. This is because the material comprising 
the Plasmodium spontaneously forms networks spanning the nutrient sources. 
Even when the plasmodium is arranged initially as a solid sheet of material, the 
sheet is soon transformed into a network structure by competitive flux of mate- 
rial within the sheet [39]. It is physically impractical to force a freely foraging 
Plasmodium to conform to a TSP network structure during its nutrient foraging, 
as shown in Fig. [T] 

Nevertheless, the material computation embodied within Physarum presents 
interesting possibilities towards generating novel spatially represented methods 
of unconventional computation. The motivation for this is threefold. Firstly, 



many natural systems exhibit properties which are not found in classical comput- 
ing devices, such as being composed of simple and plentiful components, having 
redundant parts (i.e. not being dependent on highly complex units), and showing 
resilient, or fault tolerant, behaviour. Secondly, unconventional computation is 
often observed in systems which show emergent behaviour, i.e. behaviour which 
emerges from the interactions between simple component parts, and which — 
critically — cannot be described in terms of the lower level component interac- 
tions. Emergent behaviour is characterised by systems with many simple, local 
interactions and which display self-organisation — the spontaneous appearance 
of complexity or order from low-level interactions. The study of these properties 
is useful not only from a computational perspective, but also from a biological 
viewpoint - since much of the complexity in living systems appears to be built 
upon these principles. The third reason for interest in unconventional compu- 
tation is because, for a number of applications at least, utilising the natural 
properties of physical systems is a much more efficient means of computation. 
Natural computation can take advantage of parallel propagation of information 
through a medium (for example in the chemical approximation of Voronoi dia- 
grams ([T^), or the parallel exploration of potential path choices in path prob- 
lems using microfluidic gas discharge plasmas ([M]). Stepney has suggested that 
this material computation, in utilising the natural physical properties of materi- 
als (not necessarily living systems) , may afford rich computational opportunities 

An in-silico attempt at reproducing the pattern formation and adaptation 
behaviour of Physarum using multi-agent transport networks was introduced 
in |28| and its approximation of TSP paths for very simple data layouts was 
noted. The general pattern formation abilities of the approach was characterised 
in |26| . Attempts at encouraging these multi-agent transport networks to con- 
form to TSP-like requirements (degree of connectivity 2, no crossed paths) by 
dynamically adjusting the concentration of simulated nutrient attractants, using 
a feedback mechanism based on the current configuration of the network, were 
presented in |:29 . This resulted in extremely complex transitions of network dy- 
namics and partial success in constructing TSP tours. In the approach outlined in 
this paper we attempt a simpler approach which utilises a larger aggregate mass 
of the same multi-agent collective which behaves as a morphologically adaptive 
cohesive 'blob' of virtual material. 

3 A Material Approach to TSP by a Shrinking Blob 

In the shrinking blob method we use a piece, or 'blob' of a virtual plasmodium 
material to approximate the TSP. The material is composed of thousands of 
simple mobile multi-agent particles interacting together in a 2D diffusive lattice. 
Each particle senses the concentration of a generic 'chemoattractant' diffusing 
within the lattice and each agent also deposits the same substance within the 
lattice upon successful forward movement. The multi-agent population collec- 
tively exhibits emergent properties of cohesion and shape minimisation as a 
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Fig. 1. Foraging plasmodium of Physarum does not approximate the TSP in 
both unconstrained and constrained environments, (a) Physarum plasmodia are 
inoculated at oat flakes on non-nutrient agar, (b) individual plasmodia extend 
from oat flakes and fuse, (c-d) the plasmodium continues to forage and the 
shape of the TSP is not represented, (e) foraging of plasmodium is constrained 
by placing a ring of saline soaked thread (Ig NaCl / lOOg water) at the periphery 
of the arena, (f ) as the salt diffuses into the agar the shape of the plasmodium is 
confined, (g) the pattern of the plasmodium at 36h is confined but is disconnected 
from outer nodes (note empty tube remnants) and does not approximate the 
TSP. 



results of the low-level particle interactions. The pattern formation and network 
adaptation properties of small populations of the material were discussed in [2 6) 
and were found to reproduce a wide range of Turing-type reaction-diffusion pat- 
terning. In this paper we use a relatively large population of particles which 
collectively behaves as a sheet of deformable virtual material. A full description 
of the virtual material method is given in the Appendix and an overview of the 
method follows. 

3.1 Shrinkage Process 

We initialise a sheet of the virtual material around a set of data points cor- 
responding to TSP city nodes (Fig. [2^). Chemoattractant is projected into the 
diffusive lattice at node locations, however, projection is reduced at regions which 
are covered by the blob sheet. The initial shape of the sheet corresponds to the 
Convex Hull of the data points. We then shrink the material by systematically 
removing some of its constituent particle components. The city nodes act as 
attractants to the material, effectively 'snagging' the material at the locations 
of uncovered nodes and affecting its subsequent morphological adaptation. As 
the material continues to shrink its innate minimising properties conform to the 
locations of the city nodes and the area occupied by the material is reduced, be- 
coming a concave area covering the nodes (Fig. [2jD-e) . The shrinkage is stopped 
when all of the nodes are partially uncovered by the sheet (Fig. ^) . The reader 
is encouraged to view the supplementary video recordings of the shrinkage pro- 
cess at http://uncomp.uwe.ac.uk/jeff/material_tsp.htm. The adaptation 
of the blob to the data stimuli is not entirely smooth, the video recordings show 
that the blob sheet adapts to the changing stimuli as data nodes are temporarily 
uncovered and re-covered by the blob. When the shrinkage is halted the area of 
the sheet corresponds to the area enclosed by a tour of the Euclidean Travelling 
Salesman Problem. The exact tour formed by the blob can be elucidated by 
tracking along the perimeter of the blob, adding a city to the tour list when it 
is first encountered. The tour is complete when the start city is re-encountered. 
The approach is simple, making use of the innate adaptive emergent properties 
of the material. Despite being completely unguided and containing no heuristic 
optimisation strategies the approach yields efficient tours. The separate stages 
of the approach will now be described in detail. 

3.2 Halting the Computation 

It is important to halt the shrinkage of the blob at the right time. If the shrinking 
is stopped too early an incomplete tour will be formed (i.e. only a partial subset 
of the nodes will be included in the tour if not all of the nodes are uncovered). 
Unlike guided heuristic methods a set of candidate tours is not initially formed 
and subsequently modified. Only a single tour is formed and the shrinking blob 
approach is akin to the 'instance machines' (as opposed to universal machines) 
proposed by Zauner and Conrad [15] ■ To automatically halt the computation we 
use a so-called 'traffic light' system. At the start of the method the sheet covers 



(a) t=60 



(b) t=60 



(c) t=520 




(d) t=3360 (e) t=6360 (f) t=7102 



Fig. 2. Visualisation of the shrinking blob method, (a) sheet of virtual material 
initialised within the confines of the convex hull (grey polygon) of a set of points. 
Node positions are indicated by circles. Outer partially uncovered nodes are light 
grey, inner nodes covered by the sheet are in dark grey, (b-e) sheet morphology 
during shrinkage at time 60, 520, 3360, 6360 respectively, (f) shrinkage is stopped 
automatically when all nodes are partially uncovered at time 7102. 



the entire set of nodes. Only the outer nodes are partially covered by the blob. 
To measure whether a node is covered by the sheet we assess the number of 
particles in a 5x5 window around each node. If the number of particles is < 15 
then the node is classified as uncovered and the node indicator is set to green. 
Otherwise the node is classified as covered and the node indicator is set to red. 
At each scheduler step the indicators of all nodes are checked. When all nodes 
are set to green, all nodes underneath the blob are partially uncovered and the 
shrinkage is stopped. 

3.3 Reading the Result of the Computation 

To trace the path of cities in the tour discovered by the blob a manual process 
is used. The collection of partially uncovered nodes and blob shape may be 
interpreted as an island shape with the nodes representing cities on the coastline 
of the island (Fig. [s]). We begin by selecting the city at the top of the arena. 
If more than one city is at this y location the left-most city at this y location 
is selected. This city is the start city of the tour and is added to the tour list 
T. Moving in a clockwise direction we trace the perimeter of the blob (walking 
around the shore of the island . . . ). Each time we encounter a city, it is added 
to T. If a city is subsequently re-encountered (as in the case of narrow peninsula 
structures as described below) it is ignored. When the path reaches the starting 
city the tour is complete and the list in T represents the tour of the TSP found 
by the shrinking blob. 

Some special cases in the tracking process must be noted in the case where a 
city lies on a narrow 'peninsula' of the blob as indicated in Fig. [4j In Fig. [4^ the 
city nearest position x in the path lies close to one side of a narrow peninsula. 
However the side at which the city is located can be deduced by a small convex 
bulge on the left side of the blob. In this case the city is not added until it is 
encountered on the left side of the peninsula. In the case of Fig. ^p, however, 
the city at x is located exactly in the middle of a peninsula and its closest side 
cannot be discerned. In this instance two interpretations are possible and the 
subsequent differences in possible tour paths are indicated by the dotted lines 
in Fig. [4j3, i) and ii). In interpretation i) the city is added to T immediately 
and in ii) it is not added until it is encountered on its opposite side. If this 
situation occurs during the tracking process we add the city to T when it is first 
encountered. 

4 Results 

We assessed the shrinking blob method by generating 20 datasets, each consisting 
of 20 randomly generated nodes within a circular arena in a 200x200 lattice. 
To aid the manual tracking process we added the condition that points must 
have a separation distance of at least 25 pixels. For each run a population of 
particles was generated and initialised within the confines of the convex hull 
(algorithmically generated) of the point set. Any particles migrating out of the 




Fig. 3. Reading the TSP tour formed by the shrinking blob by perimeter track- 
ing, (a) Tracking is initiahsed at the top most node. Perimeter of blob is traced 
in a clockwise direction. Each time a node is encountered for the first time it is 
added to the tour. The tour is completed when the start node is re-encountered. 




(b) 

Fig. 4. Special cases of when nodes are located on a narrow peninsula, close, 
or equidistant from either side of the 'land', (a) The node at 'x' is close to the 
middle of a narrow portion of the blob. The slight convex bulge in the blob 
indicates that it is closest to the left side and the node is not added to the tour 
until it is encountered on the left side, (b) The node at 'x' is directly in the 
middle of a narrow portion of the blob. Two potential tours are possible, shown 
in i) and ii) with their respective tours as dotted lines. If this case occurs, the 
node is added to the tour the first time it is encountered, as in i). 



convex hull area were reraoved. As the shrinkage process started the cohesion of 
the blob emerged and, as shrinkage progressed, the blob adapted to the shape of 
the city nodes. Six experimental runs were performed on each dataset and the 
resulting blob shape was recorded and tracked by the manual tracking process to 
reveal the tour. The best, worst and mean performance over the 6 runs for each 
20 datasets was recorded and these results were aggregated over the 20 datasets 
and shown in Fig. [5] Results of the shrinking blob method are expressed as a 
fraction of the shortest exact tour found by the Concorde TSP solver [TU] . 



Dataset 


Run 1 


Run 2 


Run 3 


Run 4 


Run 5 


Run 6 


min 


max 


mean 


Exact Tour 


Best 


Worst 


Mean 


Length 


Fraction 


Fraction 


Fraction 


1 


681 


683 


688 


683 


683 


683 


681 


688 


683 50 


649 


1.049 


1.060 


1.053 


2 


688 


711 


712 


735 


687 


742 


687 


742 


712.50 


636 


1.080 


1.167 


1.120 


3 


664 


664 


664 


664 


664 


664 


664 


664 


664.00 


652 


1.018 


1.018 


1.018 


4 


672 


672 


672 


672 


672 


672 


672 


672 


672.00 


665 


1.011 


1.011 


1.011 


5 


754 


717 


717 


743 


790 


717 


717 


790 


739.67 


669 


1.072 


1.181 


1.106 


6 


694 


694 


694 


670 


694 


694 


670 


694 


690.00 


659 


1.017 


1.053 


1.047 


7 


684 


762 


762 


762 


762 


684 


684 


762 


736.00 


635 


1.077 


1.200 


1.159 


8 


670 


723 


683 


723 


723 


729 


670 


729 


708.50 


667 


1.004 


1.093 


1.062 


9 


636 


636 


646 


636 


636 


636 


636 


646 


637.67 


622 


1.023 


1.039 


1.025 


10 


725 


725 


682 


682 


682 


682 


682 


725 


696.33 


674 


1.012 


1.076 


1.033 


11 


747 


748 


761 


742 


742 


742 


742 


761 


747.00 


676 


1.098 


1.126 


1.105 


12 


716 


713 


713 


720 


716 


711 


711 


720 


714.83 


684 


1.039 


1.053 


1.045 


13 


719 


682 


699 


682 


689 


699 


682 


719 


695.00 


663 


1.029 


1.084 


1.048 


14 


692 


713 


699 


701 


701 


711 


692 


713 


702.83 


623 


1.111 


1.144 


1.128 


15 


717 


692 


692 


692 


685 


685 


585 


717 


693.83 


661 


1.036 


1.085 


1.050 


16 


673 


673 


673 


673 


673 


673 


673 


673 


673.00 


670 


1.004 


1.004 


1.004 


17 


724 


714 


730 


711 


730 


711 


711 


730 


720.00 


669 


1.063 


1.091 


1.076 


18 


696 


596 


672 


571 


671 


672 


671 


696 


679.67 


630 


1.065 


1.105 


1.079 


19 


683 


703 


703 


703 


703 


703 


683 


703 


699.67 


675 


1.012 


1.041 


1.037 


20 


711 


691 


711 


711 


711 


691 


691 


711 


704.33 


644 


1.073 


1.104 


1.094 
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Fig. 5. Results of shrinking blob method over 6 runs on each of 20 randomly 
generated datasets of 20 points compared to exact results from the Concorde 
TSP solver. 



4.1 Construction of Tour by Concave Insertion Process 

Although the final tour list is read off by tracking the perimeter of the shrunken 
blob, the construction of the tour actually occurs by an insertion process as the 
blob shrinks. The blob is initially patterned with the shape of the convex hull. 
This is only a partial tour, since only the peripheral nodes which are part of the 
Convex Hull are included. By recording the stages by which nodes are uncovered 
and added during the shrinkage process, the method of construction can be 
elucidated. Fig. |8] shows the visual deformation of the Convex Hull structure as 
the blob shrinks and new city nodes are added to the list. Note that the blob 
shrinks simultaneously from all directions and the order of insertion is related to 
both the proximity of the point from the periphery of the blob and the distance 
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Fig. 6. Examples of good performance by the shrinking blob method. (a,c,e,g) 
Final blob shape with TSP tour overlaid, tour length and fraction of exact tour 
in parentheses, (b,d,f,h) Minimum exact tour found by the Concorde TSP solver. 
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Fig. 7. Examples where shrinking blob does not approximate TSP well. (a,c,e,g) 
Final blob shape with TSP tour overlaid, tour length and fraction of exact tour 
in parentheses, (b,d,f,h) Minimum exact tour found by the Concorde TSP solver. 



between two outer stimuli at the current periphery of the blob where a concavity 
forms (discussed further in Section [43] ) . The actual order of insertion of cities 
in this example is given in Fig. [9] 




(e) (f) (g) (h) 




(i) 0) (k) (1) 

Fig. 8. Construction of TSP tour by shrinking blob includes the transformation 
between Convex Hull and Concave Hull, (a) Initial Convex Hull of dataset 16 
(shown by points connected by path) is deformed to a concave shape by the 
shrinking blob. The stepwise construction of the tour is indicated by adding a 
circled point as each new city is discovered, (b-k) As the blob continues to shrink 
new points are included (circled) further reducing the area of the Concave Hull, 
(1) when shrinkage stops the set of encompassed points is a tour of the TSP. 



As the blob shrinks, concavities form in the periphery of the blob which move 
inwards to the centre of the blob shape. The concave deformation is a transfor- 
mation of the Convex Hull (CH) into a Concave Hull (OH). The Concave Hull, 
the area occupied by — or the 'shape' of — a set of points is not as simple 



M 






T 


n 


K 


G 




q 


iV 






IVl 






1 




K 


G 




Q 
O 


i\ 




R 


M 


I 




T 


D 


K 


G 




s 


N Q 




R 


M 


/ 




T 


D 


K 


G 


J 


s 


N Q 




R 


M 


I 


H 


T 


D 


K 


G 


J 


s 


N Q 




R 


M 


I 


H 


T 


D 


K 


C G 


J 


s 


N Q 




R 


M 


I 


H 


T 


D 


K 


C G 


J B 


s 


N Q 




R 


M 


I 


H 


T E 


D 


K 


C G 


J B 


s 


N Q 




R 


M 


I 


H 


T E 


D 


K 


C G 


J B 


s 


N Q 


F 


R 


M 


I 


H 


T E 


D 


K 


C G 


J B 


s 


N Q 


F 


L R 


M 


I 


H 


T E 


D 


K P 


C G 


J B 


s 


N Q 


F 


L R 


M 


I 


H 


T E 


D 


K P C G 


J B 


s 


N Q 


F 


AL R 


M 


I 


H 


T E 


D 


K P C G 


J B 


s 


N Q 


F O A L R 



Fig. 9. Gradual construction of tour by city insertion during shrinkage process 
in dataset 16, as visualised in Fig. [8] Top row shows initial configuration of blob 
as convex hull. Each row inserts a new city (bold) into the tour as indicated in 
Fig-i 



to define as its convex hull. It is commonly used in Geographical Information 
Systems (GIS) as the minimum region (or footprint Jj)]) occupied by a set of 
points, which cannot, in some cases, be represented correctly by the convex hull 
[15] . The Concave Hull is related to the structures known as a-shapes pT|. The 
a-shape of a set of points, P, is an intersection of the complement of all closed 
discs of radius 1/a that includes no points of P. An a-shape is a convex hull 
when a — oo. When decreasing a, the shapes may shrink, develop holes and 
become disconnected, collapsing to P when a — )■ 0. A concave hull is non-convex 
polygon representing area occupied by P and the concave hull is a connected 
a-shape without holes. In contrast to a-shapes, the blob (more specifically, the 
set of points which it covers) does not become disconnected as it shrinks. As the 
blob adapts its morphology from Convex Hull to TSP is demonstrates increased 
concavity with decreased area. Although the shrinkage process is automatically 
stopped when a TSP tour is formed, the process could indeed continue past the 
TSP. If shrinkage continues then the blob (now adopting a network shape) will 
approximate the Steiner minimum tree (SMT), the minimum path between all 
nodes. As demonstrated in [57] the additional Steiner nodes in the SMT may be 
removed by increasing the attractant projection from the data nodes. The ma- 
terial adapts to the increased attractant concentration by removing the Steiner 
nodes to approximate the Minimum Spanning Tree (MST). 



4.2 Blob TSP Tour as a Waypoint in the Transition Prom Convex 
Hull to Spanning Tree 



The insertion process of adding nodes to the Convex Hull reveals an orderly 
transition to the TSP which continues after further shrinkage, leading to the 
following finding. 



Finding 1 The evolution of the blob shape by morphological adaptation is a 
transition from CH to OH to TSP to MST to SMT. 



We do not explicitly include a-shapes in this transition since a-shapes can 
include holes and disconnected structures, which do not form in a defect-free 
shrinking blob. This transition is based on increasing concavity and decreasing 
area, and encompasses the a blob TSP tour bTSP as part of the hierarchy. Note 
that the blob tour bTSP is only one instance of the set of possible TSP tours 
TSP and is not guaranteed to be the minimal tour. The blob TSP tour is only 
a transient structure — a waypoint — in the natural shrinkage process (we halt 
the computation at this point merely because we are interested for the purposes 
of this report). 

It is known from Toussaint that there is a hierarchy of proximity graphs 
(graphs where edges between points are linked depending on measures of neigh- 
bourhood and closeness) [47]. Each member of the hierarchy adds edges and 
subsumes the edges of lower stages in the hierarchy and some common graphs 
(see Fig. 10 i-e) include the Delaunay triangulation DTN to Gabriel Graph GG 
to Relative Neighbourhood Graph RNG to Minimum Spanning Tree MST. 
Also shown is the shortest possible tree between all nodes formed by adding 
extra Steiner nodes (Fig. \lO^ ) ■ It was found in [2] that Physarum approximates 
the Toussaint hierarchy of proximity graphs as it constructs transport networks 
during its foraging and it was demonstrated in |27| that multi-agent transport 
networks mimicking the behaviour of Physarum also minimise these proximity 
graphs by following this hierarchy in its downwards direction. From a biological 
perspective traversing the Toussaint hierarchy suggests a mechanism by which 
Physarum can exploit the trade-off between foraging efficiency (many network 
links) and transport efficiency (fewer but fault tolerant transport links). This 
mechanism, may also be present in terms of maximising foraging area searched 
(exploration) and minimising area for efficient transport (exploitation), as sug- 
gested in . We suggest that the hierarchy we observed in the shrinking blob 
from CH to OH to TSP to MST to SMT may encompass such an area-based 
exploration-exploitation mechanism (Fig. [lO^-l) . It is notable that there is some 
overlap between the Toussaint hierarchy and the shrinking blob hierarchy where 
deepening concavities in the blob hierarchy appear to correspond to the deletion 
of outer edges in the Toussaint hierarchy, suggesting that there may be some for- 
mal relationship between the two. This possible relationship may suggest further 
studies. 
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Fig. 10. Comparison of the Toussaint Hierarchy (top row) and Shrinking Blob 
Hierarchy (bottom row), (a) Initial source nodes, (b) Delaunay Triangulation 
(DTN), (c) Gabriel Graph (GG), (d) Relative Neighbourhood Graph (RNG), (e) 
Minimum Spanning Tree (MST), (f) Steiner Minimum Tree (SMT), (g) Initial 
blob is patterned as a Convex Hull (CH), (h-i) as the blob shrinks it adopts the 
Concave Hull, (j) after uncovering the last node a TSP tour is formed, (k) blob 
can be forced to adopt MST by increasing node concentration, (1) the 'natural' 
end point of blob shrinkage is the SMT. 



4.3 Variations in Performance of the Shrinking Blob Method 

The results of the shrinking blob method show variations in performance from 
very good approximations of close-to-minimum tours (Fig. |6| to less favourable 
fractions of the minimal tour (Fig. [7|. What is the reason for the disparity 
in performance on these datasets? If we examine the tour paths we can glean 
some clues as to the difference in performance. In the 'good' results examples 
the major concave regions of the tour formed by the blob closely match the 
concavities in the exact computed TSP tour (e.g. Fig. |6^ and b). However in the 
'poor' approximation results we can see that the major concave regions of the 
blob tour do not match the major concavities in the respective exact computed 
tours (e.g. Fig. [T^ and b). Given that these concave regions are formed from the 
deformation of the initial Convex Hull we can see that the concavities in the 
blob tour appear to be formed, and deepened where there are larger distances 
between the cities on the initial Convex Hull. 

To explore the role of distance on concavity formation we patterned a blob 
into a square shape by placing regularly placed stimuli around the border of a 



square (Fig. Ill, stimuli positions, 20 pixels apart, indicated by crosses). When 
shrinkage of the blob was initiated there is no difference between the stimuli 
distances. All regions between stimuli initially show small concavities (the 'per- 
forations' in Fig. [1TJ3) until one gradually predominates and extends inwards. 
Also of note is the fact that when one concave region predominates, the other 
concavities shrink (Fig. [Tl|;-e) . The position of the initial dominating concavity 



is different in each run (presumably due to stochastic influences on the collective 
material properties of the blob) and this may explain the small differences in 
performance on separate runs using the same dataset. 

When there is a larger gap between stimulus points the predominating con- 
cavity forms more quickly and is larger. This is shown in Fig. |llf -j which has a 
gap of only 30 pixels between neighbouring stimulus points on the right side of 
the square and in Fig. |ll[ t-o which has a gap of 60 pixels between neighbouring 
stimulus points. The shorter distance between points in (f-j) generates more ten- 
sion in the sheet, prevent its deformation. The larger distance between stimuli 
in (k-o) results in less tension in the sheet at this region and the sheet deforms 
to generate the concavity. 

When there are multiple instances of large distances between stimuli there 
is competition between the concave regions and the larger region predominates. 
This is demonstrated in Fig. [TTJa-t which has a distance gap of 40 pixels on the 
left side of the square and 60 pixels on the right side. Although two concave 
regions are formed, the larger deepens whilst the smaller concave region actually 
shrinks as the blob adapts its shape. 

The synthetic examples illustrate the influence of city distance on concavity 
formation and evolution and these effects are more complex when irregular ar- 
rangements of city nodes are used. This is because arrangements of cities present 
stimuli to the blob sheet when partially uncovered, acting to anchor the blob at 
these regions, and the morphological adaptation of the blob is thus dynamically 
affected by the changing spatial configuration of uncovered city nodes. In the 
examples of relatively poor approximation of the minimum tour (Fig. [t]) the ini- 
tial incorrect selection of concavities are subsequently deepened by the shrinking 
process, resulting in tours which differ significantly in both their visual shape 
and in the city order from the optimum tour. In the examples of good com- 
parative performance with the exact solver the blob tours differ only in a small 
number of nodes. 

Although outright performance is not the focus of this report, we tested the 
blob method on a randomly generated dataset of 50 nodes in a preliminary 
assessment of scalability. We found that over 6 runs the blob method found 
tours with a best of 1.064, worst of 1.099 and mean of 1.076 compared as a 
fraction of tour length to the exact minimum computed by the TSP solver. Fig. 



12 However, these results may also be subject to the variability in performance 



seen in the 20 node examples. 



5 Discussion 



We have presented a simple material-based approach to the computation of the 
Travelling Salesman Problem using a shrinking blob. The method utilises the 
emergent morphological adaptation properties of a virtual material arising from 
local interactions within a multi-agent particle system. We shrink this material 
over time and its deformation and adaptation to the projected data points yields 
a tour of the TSP. We should again emphasise that the method is notable for 
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Fig. 11. Distance between nodes affects position and speed of concavity forma- 
tion in shrinking blob. Nodes arranged in the border of square (indicated by 
crosses on leftmost images), (a-e) all nodes have identical distance of 20 pixels, 
(f-j) right side has node gap of 30 pixels, (k-0) right side node gap 60 pixels, 
(p-t) left side node gap 40 pixels, right side node gap 60 pixels. 
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Fig. 12. The performance of the shrinking blob method on a larger dataset. (a-e) 
evolution of the shrinking blob on 50 node dataset, (f) tour formed by shrinking 
blob, (g) exact minimum tour by Concorde solver. 

its simplicity and novelty rather than its performance. Indeed the performance, 
when compared to exact TSP solvers or leading heuristic methods, compares 
relatively unfavourably in terms of absolute tour distance. The method does, 
however, contain a number of properties that are intriguing. Firstly, unlike many 
other nature inspired approaches, the method is not population based and the 
blob only computes a single instance of a TSP tour. In addition, no attempt 
is made to modify or optimise the tour. The benefits of population based ap- 
proaches are that very large search spaces in combinatorial optimisation prob- 
lems can be traversed and candidate solutions can be compared in some way. 
This allows the efficient pooling of good solutions, generation of new candidate 
solutions and avoids local minima. The blob method does not contain any of 
these beneficial features. 

How then, can a shrinking blob naturally generate a good quality (albeit non 
optimal) TSP tour? The intrinsic performance of the blob is based on its (vir- 
tual) material properties which exhibit innate minimisation behaviour. Previous 
research has demonstrated the network minimisation properties of the material 
approach, reproducing phenomena seen in soap film evolution [27j and lipid nan- 
otube networks [34] . The maintenance of uniform shape during blob shrinkage 
also allows no crossing over of paths and it is known from |18| that crossing paths 
produce non-optimal tours. More specifically, the tour is constructed by insertion 
of cities into the list as concave regions are formed in the initial Convex Hull pat- 
tern and subsequently deepened. This is similar to algorithmic heuristics which. 



beginning with a Convex Hull, add cities to the list based on certain cost criteria 
PU] . [5T] . In the case of the blob, however, there is no explicit consideration of 
cost when adding cities to the tour. The mechanism of insertion selection in the 
blob (by deepening concavities) is intrinsic to its quasi-mechanical properties of 
the 'material' which are influenced by the depth of the city to the Convex Hull 
boundary and the span distance between boundary stimuli. 

Although the blob approach differs from population based nature-inspired 
heuristics it is reminiscent, in character if not in direct operation, with other 
analogue based methods. In the Elastic Net algorithm, introduced by Durbin and 
Willshaw, a circular band is initialised at the approximate centre of a pattern 
of source TSP nodes. The band is expanded iteratively whilst two forces are 
applied to points on the band which attempt to minimise distances between cities 
on the band and the overall length of the band itself fT6]. In the conceptually 
opposite approach of [B] the band is initialised on the Convex Hull and the two 
forces attempt to constrict the band whilst attracting the band towards a city. 
The tour formed by the band is then subject to a second 'non-deterministic 
improvement' algorithm to escape local minima. The main difference between 
the blob method and these 'band' approaches is that the material properties of 
the blob method are an emergent property of, and are distributed within, local 
interactions between components of the material. The computation is thus an 
embodied property of the material itself. 

Can the mechanism underlying the simple material approximation of the 
TSP in the blob approach contribute to the question of human performance on 
the TSP? MacGregor and Ormerod noted that humans produced efficient re- 
sults on the TSP [35] and this finding stimulated further research into human 
performance on the problem and possible perceptual and cognitive mechanisms. 
In their analysis of experimental findings using human subjects Ormerod and 
Chronicle noted that global perceptual influences appear to play a role in human 
approximation of the TSP |42] . MacGregor et al. suggested a model based on in- 
sertion of cities into the Convex Hull |37j . This model is similar to the shrinking 
blob mechanism except that in the blob approach the addition of cities occurs 
in parallel whereas in the MacGregor et al. model it is a sequential process. The 
blob method also exhibits another property found in optimal tours, that bound- 
ary points in the original convex hull are connected in sequence (the sequence 
may, of course, be interrupted by interior points). Other competing models to 
explain human TSP performance exist, including variants of hierarchical pyra- 
mid models [5T], [13] and the global- local model proposed by Best [TT]. Merits, 
problems, biological plausibility and the role of local vs. global perceptual pro- 
cesses of the competing models have been the subject of lively debate and an 
assessment is beyond the scope of this paper, but see the review in [35] for an 
overview. In this review MacGregor states that, despite the growing interest in 
research into the human performance on the TSP, and combinatorial optimi- 
sation problems in general: "As yet, no algorithms have been put forward to 
explain performance on the MSTP [Minimum Spanning Tree Problem] and the 
GSTP [Generalised Steiner Tree Problem] . . . " . It is notable that the shrinking 



blob method incorporates approximations of all three problems and executes 
a natural transition from global to local 'perception' using material properties 
which emerge from very simple low-level and bottom-up interactions. Whether 
this natural computation employed by the blob is of interest, or utility, to human 
combinatorial optimisation problems is, however, an open question. 

Limitations of the approach, as it stands, include its noted relatively modest 
performance and the reliance of a manual method to interpret the result of the 
blob computation. Manual interpretation of the blob tour is, of course, open to 
experimenter bias and for this reason a methodical process must be followed, as 
described in section 3.3 An automated method of tracking the perimeter and 
'reading' the result of the blob tour would, nevertheless, be of benefit. Further 
work, including a comprehensive evaluation of model parameters affecting the 
material properties of the blob may suggest methods by which the basic features 
of the shrinking blob approach may be adapted, or improved, to improve the per- 
formance in comparison with leading heuristic methods. The material properties 
and computation of the blob emerge from a population of simple multi-agent par- 
ticles and It would be satisfying if this virtual material could be implemented 
and embodied in a real physical substrate with the desired physical (for example 
visco-elastic, free energy minimisation) properties. Alternately it may be possi- 
ble to translate the material operation of the unconventional computation blob 
method into a classical algorithmic method. 
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7 Appendix: Shrinking Blob Particle Model Description 

The multi-agent particle approach to generate the behaviour of the virtual ma- 
terial blob uses a population of coupled mobile particles with very simple be- 
haviours, residing within a 2D diffusive lattice. Lattice size was 200 x 200 pixels. 
The lattice stores particle positions and the concentration of a local diffusive 
factor referred to generically as chemoattractant. Collective particle positions 
represent the global pattern of the blob. The particles act independently and 
iteration of the particle population is performed randomly to avoid any artifacts 
from sequential ordering. 



7.1 Generation of Emergent Blob Cohesion and Shape Adaptation 

The behaviour of the particles occurs in two distinct stages, the sensory stage 
and the motor stage. In the sensory stage, the particles sample their local en- 
vironment using three forward biased sensors whose angle from the forwards 



position (the sensor angle parameter, SA), and distance (sensor offset, SO) may 
be parametrically adjusted (Fig.[l3^). The offset sensors generate local coupling 
of sensory inputs and movement to generate the cohesion of the blob. The SO 
distance is measured in pixels and a minimum distance of 3 pixels is required for 
strong local coupling to occur. During the sensory stage each particle changes 
its orientation to rotate (via the parameter rotation angle, RA) towards the 
strongest local source of chemoattractant (Fig. [Tsja). After the sensory stage, 
each particle executes the motor stage and attempts to move forwards in its cur- 
rent orientation (an angle from 0-360 degrees) by a single pixel forwards. Each 
lattice site may only store a single particle and particles deposit chemoattractant 
into the lattice only in the event of a successful forwards movement. If the next 
chosen site is already occupied by another particle move is abandoned and the 
particle selects a new randomly chosen direction. 

7.2 TSP Problem Representation 

Twenty datasets were generated, each consisting of 20 randomly chosen data 
points within a circular arena. A condition was added that a minimum distance 
of 25 pixels must exist between data points. This gives a more uniform distribu- 
tion of node points, preventing clustering of node points often found in real- world 
TSP instances, for example those based on real city locations. This was done 
partly to aid visual tracking of the tour path but also because a less clustered 
node distribution is known to provide more challenging problem instances (e.g. 
in human performance on TSP [24]) since there is less likelihood of providing 
pre-existing cues to intuitive solutions (for example nearest neighbour grouping). 
Each dataset was saved to a text file. TSP city data points for each dataset were 
loaded from the text files and were represented by projection of chemoattractant 
to the diffusion lattice at locations corresponding to a 3 x 3 window centred about 
their x^y position. The projection concentration was 1.275 units per scheduler 
step. If a node was covered by a portion of the blob (i.e. if the number of agents 
in a 3 X 3 window surrounding the node was > 0) the projection was reduced 
to 0.01275 units. Suppression of projection from covered sites was necessary to 
ensure a uniform concentration within the blob at internal data points. Uncov- 
ering of the data points by the shrinking blob acted to increase concentration at 
exposed nodes, causing the blob to be anchored by the nodes. Diffusion in the 
lattice was implemented at each scheduler step and at every site in the lattice 
via a simple mean filter of kernel size 3x3. Damping of the diffusion distance, 
which limits the distance of chemoattractant gradient diffusion, was achieved by 
multiplying the mean kernel value by 0.95 per scheduler step. 

The blob was initialised by creating a population of particles and inoculating 
the population within the bounds of a Convex Hull formed by the data points, 
generate at the start of the experiment by a conventional algorithm. Why choose 
the Convex Hull as the initial blob pattern? It is possible to select any pattern 
for the blob shape, for example a circular blob. However, selecting any arbitrary 
pattern to cover the points might leave a large portion of the blob overlapping 
vacant space. This could result in shrinkage in one part of the blob taking place 
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Fig. 13. Architecture of a single component of the shrinking blob and its sensory 
algorithm, (a) Morphology showing agent position 'C and offset sensor positions 
(FL, F, FR), (b) Algorithm for particle sensory stage. 



over vacant space with shrinkage in other parts over city points and would bias 
the output result. The Convex Hull was thus chosen as a simple representation 
of the perimeter of the data points. The exact population size differed depending 
on the area of the Convex Hull but the initial blob was typically composed of 
between 10000 and 15000 particles. Particles were given random initial positions 
within these confines and random initial orientations. Any particles migrating 
out of the bounds of the Convex Hull region were deleted. Particle sensor offset 
(SO) was 7 pixels. Angle of rotation (RA) and sensor angle (SA) were both 
set to 60 degrees in all experiments. Agent forward displacement was 1 pixel 
per scheduler step and particles moving forwards successfully deposited 5 units 
of chemoattractant into the diffusion lattice. Both data projection stimuli and 
agent particle trails were represented by the same chemoattractant ensuring that 
the particles were attracted to both data stimuli and other agents' trails. The 
collective behaviour of the particle population was cohesion and morphological 
adaptation to the configuration of stimuli. 



7.3 Shrinkage Mechanism 

Adaptation of the blob size is implemented via tests at regular intervals as 
follows. If there are 1 to 10 particles in a 9 x 9 neighbourhood of a particle, 
and the particle has moved forwards successfully, the particle attempts to divide 
into two if there is a space available at a randomly selected empty location in 
the immediate 3x3 neighbourhood surrounding the particle. If there are to 80 
particles in a 9x 9 neighbourhood of a particle the particle survives, otherwise it is 
deleted. Deletion of a particle leaves a vacant space at this location which is filled 
by nearby particles, causing the blob to shrink slightly. As the process continues 
the blob shrinks and adapts to the stimuli provided by the configuration of city 
data points. The frequency at which the growth and shrinkage of the population 
is executed determines a turnover rate for the particles. The frequency of testing 
for particle division was every 5 scheduler steps and the frequency for testing for 



particle removal was every 10 scheduler steps. Since the shrinking blob method is 
only concerned with the reduction in size of the population it might be asked as 
to why there are tests for particle division at all. The particle division mechanism 
is present to ensure that the adaptation of the blob sheet is uniform across the 
sheet to prevent 'tears' or holes forming within the blob sheet, particularly at 
the start of an experiment when flux within the blob is initially established. 

7.4 Halting Mechanism 

The shrinkage of the blob was halted when all data points were partially un- 
covered by using the following calculation. At each data point the number of 
particles surrounding the point in a 5 x 5 window was sampled. If the number of 
particles in this window was < 15 the node was classified as partially uncovered. 
When all nodes are uncovered the model is halted and a greyscale representation 
of the diffusive lattice is saved to disk to read the result of the blob tour. The 
calculated tour path was saved in a text file and tour distance was calculated. 
The exact minimum tour path for each dataset was calculated for comparison 
by loading the data point configuration files into the Concorde TSP solver [TU] . 
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